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ABSTRACT 

Telescope Point Spread Function (PSF) quality is critical for realising the potential of 
cosmic weak lensing observations to constrain dark energy and test General Relativity. 
In this paper we use quantitative weak gravitational lensing measures to inform the 
jy^ \ precision of lens optical alignment, with specific reference to the Dark Energy Survey 

(DES). We compute optics spot diagrams and calculate the shear and flexion of the 
PSF as a function of position on the focal plane. For perfect optical alignment we verify 
the high quality of the DES optical design, finding a maximum PSF contribution to 
the weak lensing shear of 0.04 near the edge of the focal plane. However this can be 
increased by a factor of approximately three if the lenses are only just aligned within 
their maximum specified tolerances. We calculate the E and B-mode shear and flexion 
variance as a function of de-centre or tilt of each lens in turn. We find tilt accuracy to be 
a few times more important than de-centre, depending on the lens considered. Finally 
we consider the compound effect of de-centre and tilt of multiple lenses simultaneously, 
by sampling from a plausible range of values of each parameter. We find that the 
compound effect can be around twice as detrimental as when considering any one lens 
alone. Furthermore, this combined effect changes the conclusions about which lens is 
most important to align accurately. For DES, the tilt of the first two lenses is the most 
important. 

Key words: telescopes, surveys, gravitational lensing 



> 
o 

<N 
CO 



o 



1 INTRODUCTION (HSC) survejQ and the Dark Energy Survey: DE!©, tele- 

scopes under construction (the Large Synoptic Survey Tele- 
Weak lensing cosmic shear has great potential to be one 
of the most po werful tools available t o uncover the nature 
of da rk energy (|Albrecht et al.l 120061 ; IPeacock &; Schneider! 
2006). A number of planned and forthcoming surveys 

plan to use this probe of cosmology, including imminent i http://www.naoj.org/Projects/HSC/HSCProject.html 

surveys (Kilo-Degree Survey: KIDS, Hyper Suprime-Cam 2 http://www.darkenergysurvey.org 
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Figure 1. Diagram of the optical layout for DECam, shown in the 
orientation it will take when at the prime focus on La Blanco. Two 
light rays are coming from the primary mirror below the optical 
corrector. The blue light ray is an on- axis beam and the green 
ray is an off-axis beam. In our study, the lenses are decentred in 
the y direction and tilted around the z direction. 

scope: LSST0), and future space telescopes (EuclicQ and 

wfirstB). 

In the case of DES, weak lensing is one of four indepen- 
dent methods which will be used to determine the dark en- 
ergy equation of state parameter, w, to a precision of better 
than 5%. The three other methods are galaxy cluster sur- 
veys, galaxy angular clustering and super novae lig ht curves. 
Details have been described in lAnnis et al.l (| 20041 ). 

In weak lensing, galaxy shapes are distorted by the cur- 
vature of intervening space-time caused by matter in the 
Universe. As the majority of matter is dark and therefore is 
difficult to see by other means, lensing provides a useful tool. 
The gravitational lensing effect can usually be described by 
a slight squashing of galaxy images, called shear, and a slight 



3 http://www.lsst.org 

4 http://sci.esa.int/euclid 

5 http://wfirst.gsfc.nasa.gov 



bending, called flexion. The vast majority of effects are ex- 
tremely small; for example an intrinsically circular galaxy 
would typically be sheared into an ellipse with a major to 
minor axis ratio of about 1.01. 

Weak lensing directly probes the gravitational potential 
along a line of sight. Unfortunately, the atmosphere and tele- 
scope imaging also distort galaxy images, and this is usually 
a much bigger effect than the gravitational lensing effect we 
are trying to measure. The most important effect of the at- 
mosphere and telescope is well described by a convolution of 
the image with the point spread function (PSF). The PSF 
can be measured from images of point sources (stars); if the 
PSF is well known then it can in principle be removed, al- 
lowing a noisy estimate of the galaxy shear and flexion to 
be recovered. However, in practice it often leaks into the 
lensing measurement s (e.g. due to model bias or noise bias 
IVoigt fe Bridlell2010l : ?). 

A key concern in developing the Dark Energy Survey 
programme, then, is ensuring that the Dark Energy Sur- 
vey camera, DECam, has an optical system which does not 
add substantial image distortions. DECam has been con- 
structed, with careful optical alignment carried out. Instal- 
lation is now taking place during the first half of 2012, and 
the survey is due to start in Autumn 2012. The large lens 
size and weight make building work challenging due to the 
tight tolerances on the positioning. These tolerances have 
been designed to keep the instrumentation distortions to a 
minimum, so that their contributions to the weak lensing 
systematics is minimal, but even so, the expected ellipticity 
of the PSF could be larger than the lensing signal, so careful 
modelling of the optical distortions will be needed to recover 
cosmological information. 

In order to ensure that aspects of the alignment have 
been prioritised in proportion to the influence they will have 
on weak lensing systematics, we have engaged in the studies 
presented in this paper. We examine what impact errors in 
lens alignments have on the camera contribution to the lens- 
ing signal. The work shown here is specific to the DECam 
optical design, but also acts as a case study for development 
of future lensing-optimised optical telescopes. This work fol- 
lows on from e a rlier w ork on PSF requirements from lensing 
by I Kent et alj (|2006h , who computed PSF distortions for 
DES for a variety of different contributions from the optics 
and compared the result with those from other telescopes. 

This paper is organised as follows. In section [2] we re- 
view the design of DECam. We discuss the shear and flexion 
measures used to evaluate the optical performance in sec- 
tion [3] We present our results in section [4J and discuss the 
implications for lensing measurements in section [5] 



2 OPTICAL DESIGN OF DECAM 

DECam is a new wide-field prime focus camera. The optics 
consist of five lenses, ranging in diameter from 980mm to 
542mm (|Antonik et al]l2009l : [Flaudierll2010h . with six filters 
between the third and fourth lenses covering wavelengths of 
350 to llOOnm, the layout of which is shown in flgure[TJ Lens 
5 acts as the window to the cryostat that holds the CCDs; 
this cooling red uces CCD noise. Lenses 1 to 4 are held in 
nickel-iron cells (|Doel et al . 2008) glued into position using 
a silicone rubber solution. Lens 5 is held in its cell during 
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Figure 2. The camera point spread function at the centre of the focal plane (top left), and at the edge of the focal plane (bottom left). 
Right panels show the corresponding PSFs combined with a 0.7" FWHM Gaussian, to describe the PSF including atmospheric seeing. 
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Table 


1. De-centre 


and alignment 


tolerances of the lenses for 



assembly and stability during operation, from Doel et al. 2008. 



Centre Thickness Edge Thickness Diameter Weight 
Lens (mm) (mm) (mm) (kg) 



CI 


112.1 


74.34 


980 


172.36 


C2 


51.285 


148.11 


690 


87.20 


C3 


75.1 


38.27 


652 


42.62 


C4 


101.68 


52.48 


604 


49.69 


C5 


54.68 


36.19 


542 


24.37 



Table 2. Dimensions of the DECam lenses, from Doel et al. 2008. 



operation by the vacuum under which it is placed. These 
cells are attached to the barrel which is supported in the 
primary cage above the primary mirror. 

The construction of the optical corrector includes two 



main alignment stages: (i) the alignment of the lenses into 
their cells and (ii) the alignment of the lens-cells into their 
holding barrel, which is held at prime focus. In order for 
precisely focused images to be created, each part of the op- 
tical corrector must be aligned to within a series of carefully 
denned tolerances including (i) manufacturing tolerance; (ii) 
"assembly tolerance" to which the element must be aligned 
with respect to the nominal optical axis of the camera, when 
the optical corrector is being built; and (iii) the acceptable 
"dynamic tolerance" amount by which the lenses can move 
as the optical corrector is moved across the sky. The optical 
corrector is designed to minimise movement once assembled; 
however, such motions cannot be entirely suppressed and so 
during operation the optical alignment will change slightly 
depending on the position of the telescope relative to the 
zenith. Tolerances to which the lenses must be aligned are 
given in Table [T] These tolerance are challenging to meet due 
to the size and weight of the lenses, dimensions of which are 
given in Table [2] 

This paper focuses on the effect of non-perfect align- 
ment of the lenses which occurs during the lens-cell assem- 
bly stage, when the lenses are glued into their holding cells. 
While lenses were not glued until they are within tolerance, 
perfect alignment is not possible within any realistic time 
frame. It has therefore been useful to discover which lenses 
affect the weak lensing data the most, to ensure higher pre- 
cision of alignment was attempted for these lenses. In the 
following sections we explore which lenses need the most 
attention to increase the science output of the camera. 
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Figure 3. Whisker plots for the shear, 1- and 3- flexion, distributed over the full focal plane. Bottom left corner shows the scale of the 
distortion; flexion is in units of arcsec - 1 . 



3 SHEAR AND FLEXION 



In this section we describe the requirements and goals we 
set for the telescope and atmospheric point spread func- 
tions (PSF). We assume throughout that the atmosphere 
and CCDs contribute a circular Gaussian component of Full- 
Width Half Maximum (FWHM) 0.7 arcseconds, represent- 
ing good seeing conditions at Cerro Tololo where the median 
with the Blanco t elescope together with the Mosaic imager 
is 0.9 arcseconds (|Annis et al.ll2004l ). We convolve this with 
the optics point spread function to obtain the full PSF which 
we consider. 

To perform a weak lensing analysis, a galaxy image 
is corrected to remove the impact of the PSF as well as 
possible. However, several effects can cause an imperfect 
PSF rem oval, for example insufficient knowledge of the 
PSF fe.g. IPaulin-Henriksson et~a l. 2008) and a wro ngly as- 
sumed model for the galaxy (|Voigt fe Bridi3l2Q10h . These 
PSF residuals will propagate into the galaxy shear es- 
timates and contaminat e the cosmological analysis (e.g. 
lAmara fe Refregierll2007l ). The amount by which the resid- 
uals propagate depends on the galaxy shear measurement 
method used, as demon strated in a series of weak lensing 
data analysis chall enges ([Hevmans et al.ll2006l ; lMassev et al.l 
120071 : [Bridle 2010). Therefore in this paper we mainly con- 
sider the raw PSF ellipticity shear and flexion measures, and 
only consider later the probable fraction of this which might 
leak into weak lensing estimates. 



We generate PSF convolution contributions for the at- 
mosphere and the instrument separately. A simple circular 
Gaussian is used for the atmosphere. The instrument PSF 
depends on the optical layout, which we input into the opti- 



cal modelling program ZEMA>s0- This program allows sim- 
ulated light rays to be traced through an optical system and 
produces a resulting PSF at a given position in the focal 
plane. For all the following results the central ray of the i 
band filter, at a wavelength of 743. 9nm, was used. We de- 
fer discussion of the impact of the wavelength depend ence 
of the DES PSF (|Cvpriano et al.ll2010l : Ivbigt et aUl2012h to 
future work. At all times the focal plane remained in the 
same position at perfect alignment. 

Example images of the generated PSFs are shown in 
Figure[2j before and after they are combined with a Gaussian 
atmosphere, at positions 0,0 and -0.54,0.46 degrees from the 
centre of the field. Moments were then used to convert the 
PSFs into shear and flexion measures using first moments 



f I(x,y)xdxdy 

x — ,J - 1 . — 

J I(x,y)dxdy 

_ j I(x,y)ydxdy 

J I(x,y)dxdy 

quadrupole moments 

j /(x, y) (x — x) 2 dx dy 



Q: 



Q x 



j I(x,y) dxdy 
J I(x,y)(x-x)(y- y) dx dy 
j I(x,y) dxdy 

J i(x,y) (y -y) 2 dxdy 



(1) 
(2) 

(3) 
(4) 
(5) 



yu J I{x,y) dx dy 

and the further moments given by lOkura et al.l (|2008l ) re- 
quired for 1- and 3-flexion estimates. The ellipticity can then 
be found for a noise-free image using 

6 http://www.zemax.com 
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Figure 4. Convergence maps for the shear (left panels), 1-flexion (middle panels) and 3-flexion (right panels) for a perfectly aligned 
optical system. Upper panels show the E modes and lower panels show B modes. Flexion is in units of arcsec - 1 . 



ei + ie 2 — 



Qxx Qyy X u 



Qxx ~h Qyy H~ '^'(QxxQy 



(6) 



(jBonnet fc Mellie d Il995h where we introduce the standard 
complex number notation e — a + i€2 where i 2 = — 1. An 
intrinsically circular galaxy becomes stretched into an ellipse 
with ellipticity e on application of a shear 7. We will apply 
the same notation to PSFs and use the above equation to 
calculate the 'shear' of the PSF. We discuss how this impacts 
on galaxy shear measurements in section 4. 

The correspondin g equations r e lating moments to flex- 
ions are provided by lOkura et al.l (|2008l ). equations (23)- 
(24). Similarly to shear, we use these equations to calculate 
the flexion of the PSF. 

We show example shear and flexion whisker plots in 
Fig. [3j The upper panels use the best-case alignments, and 
we can see that the shear is everywhere less than 4 per cent 
and the flexion is less than 3 per arcsec. For comparison, we 
also show the results for a "worst-case" misalignment, as de- 
fined during the optical design: that giving the largest RMS 
spot size. This has a large area with shear values greater than 
0.1, reaching 0.2 at the edges, with a significant tangential 
component which would indicate a spurious mass concentra- 
tion in the center of the field if unaccounted for. The flexion 
maps have values up to 5 per arcsec, and generally shows 
a displacement in a single direction for the 1-flexion, and a 
trefoil angled in a similar pattern to the perfect alignment 
case for the 3-flexion. 

These plots of PSF distortion as a function of position 
on the focal plane can be used to reconstruct the projected 
matter density map that would be obtained if they w e re un - 
corrected, using extensions of the iKaiser &; Squires! fl993) 
reconstruction method. This gives both real and imaginary 
density distributions, depending on the orientations of the 
distortions. This is because lensing by matter usually only 
produces particular patterns of distortions. It is frequently 
assumed that the imaginary density map from real data 
gives an indication of the amount of systematic e.g. tele- 
scope effects which have leaked into the real density map. 



It is therefore of considerable interest to examine whether 
typical telescope distortions produce similar amounts of real 
and imaginary contributions to the density map, and so we 
consider both real and imaginary densit y maps here. For 
shear the required equations are the usual IKaiser &; Squired 
(1993) results 



Ks(k) = j(k)D(k) 



(7) 



where ks is the Fourier space scaled projected density (con- 
vergence) map, 7 is the Fourier transform of the PSF shear 
map and D is the response of a delta function shear map 



D 



&2 



- 2ikik2 



(8) 



k-^ ~~\~ k<2 

We extend the Kaiser & Squires formalism to flexion and 
find 



n F (k) = F(k)Q(k) 

q _ -iki - k 2 
k-^ ~\~ k<2 

and 

K~ G (k) = G(k)R(k) 

3fcifc|) + (fcf-3fc?fc 2 ) 



R 



-i(ki 



(9) 
(10) 

(11) 
(12) 



{kl + kiy 

Fig. [4] shows the resulting convergence reconstructions for 
the best case PSF shear and flexion maps. The PSF shear 
is roughly circularly symmetric, with tangential distortions 
in the centre surrounded by radial distortions. It therefore 
makes sense that the E-mode density reconstruction (top 
left panel of Fig. |4| is positive in the centre with a sur- 
rounding negative ring. Since shear is non-locally related to 
density, the central tangential shears are also partially ac- 
counted for by positive density around the edges of the field, 
which occurs especially at the diagonal corners of the map. 
The B-mode density reconstruction from the shear map is 
an indicator of the extent to which the shear map is not cir- 
cularly symmetric. In the centre part of the field the range 
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(te shear 


<jb shear 


(te 1-flexion 


(Tb 1-flexion 


(7e 3-flexion 


(tb 3-flexion 


Perfect case 


0.0086 


0.0044 


0.0026 


2.29E-5 


0.0168 


0.0038 


Worst case 


0.0557 


0.0244 


0.0093 


0.0049 


0.0317 


0.0276 



Table 3. Standard deviation of the convergence maps for shear and flexion for perfect alignment and worst case alignment. Flexion is 
in units of arcsec -1 . 



in the B-mode density is slightly smaller than that for the 
E-mode density, as can be expected by the symmetry of the 
perfect alignment case. The non-trivial relationship between 
the size of the E-mode density and B-mode density across 
the field is particularly interesting in the context of sys- 
tematics assessments in cosmic shear, which often assume a 
similar amplitude for both modes. 

The 1-flexion translates into a density map which is an 
order of magnitude smaller than that from the shear, which 
bodes well for use of 1-flexion for analysing true density 
peaks. The density range for the 3-flexion is similar to that 
for shear, with a significant negative peak in the center of 
the field. This is expected from the radially growing trefoil 
shape of the PSF seen in the top right panel of Fig. 3. 

In order to encapsulate the convergence maps into a 
single number for comparisons between camera configura- 
tions, we use the standard deviation, cr sy s, of the conver- 
gence map. This statistic is important for statistical analy- 
ses of shear and flexion, for example cosmic shear correlation 
functions. Standard deviations are given in Table [3] for the 
perfect alignment and worst case alignments. As expected 
from the density map (Figure 4), the largest standard devi- 
ations for perfect alignment are for the E-mode shear and 
E-mode 3-flexion. For both perfect alignment and worst case 
alignment, the standard deviation for E-mode shear and B- 
mode shear are a comparable size, which suggests that sys- 
tematics checks of the B-mode are a useful diagnostic of PSF 
leakage in the E-mode signal. However, the B-modes are a 
factor of two smaller than the E-modes, even for worst case 
alignment, so the exact value of the B-mode signal cannot be 
taken as an estimate of E-mode contamination (e.g. cannot 
be simply subtracted off the E-mode to find a true E-mode). 
The values for 1-flexion are still very small even for the worst 
case alignment, while the 3-flexion becomes large for both 
E- and B- mode in the worst case. 



4 EFFECTS OF MISALIGNMENT 

In this section we examine the effect of lens misalignments 
on the image quality using two methods. We first move each 
lens in turn with the others remaining in perfect alignment. 
This allows us to observe the effect which each lens has on 
the system. Secondly, we move all the lenses by a random 
amount, whilst staying within the given tolerances, to see 
the effect of multiple misalignments. Looking at multiple 
lens displacements allows an examination of typical shear 
values which the system will exhibit. 

4.1 Relative effect of single lenses 

We first de-centre each lens in turn while keeping all other 
lenses in perfect alignment. The resultant shears and flexions 



of the PSF generated are then converted into E and B mode 
convergence maps, and we calculate standard deviations ge 
and ob- These are plotted as a function of de-centre distance 
in Fig. \5\ The x-axis is the de-centre distance in microns, 
and the lines terminate at the tolerance criteria. Note that 
the y-axis does not start at zero, and thus we note that for 
the quantities generating the largest convergence variances 
(shear and 3-flexion) the exact misalignment value of any 
one lens does not change the convergence variance by a very 
large amount. 

We see that lens CI has the biggest effect of any single 
lens which is not suprising as it is highly curved and at the 
start of the optical system. Lens C2 also has a comparable 
effect on 1- and 3-flexion, while the other lenses have an 
effect on the variance more than ten times smaller than CI. 

We next consider a tilt of the lenses in Fig. \6\ We see 
that the lens tilts are much more important for lensing mea- 
surements than the lens de-centres, especially the tilt of the 
lenses with aspheric surfaces C2 and C4. For the most impor- 
tant cosmological quantity, E-mode shear, the convergence 
variance roughly doubles as the C4 tilt goes from perfect 
alignment to the full tolerance range, while tilt of other in- 
dividual lenses has relatively little impact. The E-mode 1- 
flexion is not greatly changed from its original small value 
by the tilts, and although the B-mode first flexion increases 
by a large factor especially due to C4 tilt, the value is still 
small. C4 also has the biggest effect on the third flexion, and 
CI and C2 have considerably smaller effects. Both C3 and 
C5 have little effect when tilted or decentred. 



4.2 Multiple lens offsets 

We now wish to answer two questions: (i) Does the relative 
importance of the lenses hold when more than one lens is 
offset from the ideal position? (ii) What is a realistic size of 
distortion that will come from the completed optical cam- 
era? In order to answer these questions we sampled 100 dif- 
ferent alignments from a multi-variate Gaussian distribution 
of de-centres and tilts, where the standard deviation for each 
de-centre or tilt is equal to one third of the tolerance. 63 of 
these landed within the tolerance boundaries for the camera 
alignment, and we will consider this sample. 

The resulting convergence standard deviations are plot- 
ted as a histogram in Fig. [71 The most important quantity 
for cosmology, the shear E-mode convergence range, is shown 
in the top-left panel. It has the widest range of all the quan- 
tities, which extends by more than a factor of two beyond 
the values obtained by varying the alignment of any one 
lens. It also extends below the perfect alignment case, pre- 
sumably due to certain cancellations. A similar but slightly 
less extreme pattern exists for the B-modes from the PSF 
shears. The flexion histograms are quite tight around the 
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Figure 5. The rms real part (top) and imaginary part (bottom) of the projected density reconstructed from the PSF shape component, 
as a function of de-centre of each lens individually. The left hand panels are for shear, the center panels are 1-flexion and the right hand 
panels are 3-flexion. Flexion is in units of arcsec -1 . 
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Figure 6. As Fig. 5, but as a function of the tilt of each individual lens. 
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Figure 7. Histogram of 63 random misalignments of all lenses, showing ge and gb of the system. Left panel is shear, centre panel is 
1-flexion, and right panel is 3-flexion. Flexion is in units of arcsec - 1 . 
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Figure 8. Plots of the systematic ge shear against the tilts on CI, C2 and C4. The value in the top right hand corner is the correlation 
coefficient. Tilt is in units of mm. 



perfect alignment case, and thus place less pressure on the 
need for stringent alignments. 

Next we investigate which lens is causing the greatest 
distortions, in this multi-parameter variation. We plot the 
convergence standard deviation a sys against the tilt and de- 
centre of each lens. The most significant examples are shown 
in Fig. [8] for the misalignments which had the biggest effect 
when varied individually. We calculate the correlation coef- 
ficient between the a sys values and the tilts and decenters 
of each lens, showing the results in Table [4] By comparing 
with random samples from Gaussian uncorrelated distribu- 
tions with the same number of samples, we find that a value 
above 0.26 means that the lens gives a significant contribu- 



tion to either the shear or flexion, at 95% confidence. It is 
interesting to note that C4 tilt, which looked very impor- 
tant in figure [6] has little correlation with the a sys of the 
multiple lens off-sets. However, the tilt of lenses CI and C2 
now dominate the distortions. This is in keeping with the 
previous DECam design observation that CI and C2, when 
displaced, have the largest impact on the size of the root 
mean square spot size of the PSF. 
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C3 tilt 
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0.1744 
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Third Flexion 
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-0.0111 


C4 tilt 
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0.1363 
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First Flexion 


0.1406 
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Shear 


0.2526 


0.2218 


First Flexion 


0.1595 


0.1838 


Third Flexion 


0.2150 


0.2498 



Table 4. Table showing correlation coefficients for E-mode and 
B-mode distortions, from the covariance matrix for random mis- 
alignments, under the given tolerances. Flexion is in units of 
arcsec -1 . 



5 IMPACT ON GALAXY SHEAR 
MEASUREMENT 

In order to recover a weak lensing signal, the systematics in 
the residual shear must be controlled to allow the results to 
be statistically significant. This tolerance on the variance, 
cr 2 ys can be denned by area of the sur vey A s , the galaxy 
dens ity n g , and the median redshift Zm (|Amara &; Refregieri 
l2QQ7h where 

asys<1 ° \2x lOMeg 2 ) V35arcmin- 2 J VOgJ ( 







Residual Shear 
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jmment 


2.7xlCT 4 i.4xicr 4 


Worst alig 


;nment 


4.3xlCT 4 2.5xlCT 4 



Table 5. Expected size of remaining shear in the image after 
the telescopic distortions have been removed to the extent that 
current analysis techniques allow. 

For DES, we have A s = 50 00 deg 2 , n g ~ 10 arcmin -2 and 
z m ~ 0.7 ([Annis et al.l |2003 ) . This puts a requirement on 
DES for a sys to be less than 0.0008. While all the results 
quoted previously are significantly above this level, this has 
not yet taken into account the expected correction by the 
weak lensing pipeline. 

F rom the comparison of current methods by iBridld 
the residual shear after correction would be a factor 
of at least 120 times smaller than the PSF shear described 
in this paper. R eductions in the flexions were not calculated 
m Bridlel (|2010h . so it is not yet known whether they can be 
reduced by a similar amount. However, for shear, this fac- 
tor brings cr syss hear of the residual shears under the required 
8x10 ~ 4 for both perfect and worst case alignment, as shown 
in table [5] 



6 CONCLUSIONS 

In this paper we have explored the effect of camera optical 
alignments on quantities important for weak gravitational 
lensing, with specific reference to the Dark Energy Survey 
(DES). We have calculated the shear and flexion of the PSF 
as a function of position on the focal plane, for a realistic ray- 
tracing model of DECam including misalignments of lenses. 
In this work, we do not include flexure errors, or spacing 
tolerances which can be overcome by adjusting focus. We 
also do not consider flatness, position or tilt of the focal 
plane. However, we include lens decentres and tilts, which 
are known to be within tolerance for the actual DECam. 

For perfect optical alignment, we find a maximum PSF 
shear of 0.04 near the edge of the focal plane. However this 
can be increased by a factor of roughly three if the lenses 
are only just aligned within their maximum specified toler- 
ances. In both cases, however, once PSF correction methods 
are employed, the resulting impact on shear estimates for 
galaxies in the DES weak lensing catalogue are expected to 
be below the required threshold for measuring dark energy 
parameters. 

We have calculated the E and B-mode shear and flexion 
variance as a function of de-centre or tilt of each lens in turn. 
We found tilt accuracy to be a few times more important 
than de-centre, depending on the lens considered. Finally, 
we have considered the combined effect of de-centre and tilt 
of multiple lenses simultaneously, by sampling from the per- 
mitted range of values of each parameter. In this case we 
find that combined effect can be around twice as detrimen- 
tal as when considering any one lens alone. Furthermore, 
this combined effect changes the conclusions about which 
lens is most important to align accurately; for DES, the tilt 
-of the first two lenses is the most important. 

The results of these simulations have been used to in- 
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form alignment of the DECam lenses (i.e. on which lens 
alignments to spend the most effort), and acts as an impor- 
tant confirmation that the DECam optical tolerances lead 
to a system which is fit for extremely accurate weak lens- 
ing measurements. The final image quality of course will be 
affected by other factors such as barrel sag, optical wedges, 
surface figure errors, focal plane misalignment and material 
inhomogeneity (in roughly decreasing order of significance). 

For example, the quality of the lens polishing intro- 
duces surface figure errors. We have assumed perfect pol- 
ishing for this work. The smoothness of the final DES lenses 
have been examined using interferometry and photography 
to assess phase deviations and is well within specifications. 
Preliminary investigations suggest that typical atmosphere 
convolved PSF ellipticities are increased by up to twenty per 
cent in a way which varies in a non- symmetric way across 
the field of view. It appears to be a small but non-negligible 
fraction of the effect due to random misalignments of lenses 
within their tolerances studied in this work. 

Work is ongoing to incorporate all the image quality 
effects into the optical model along with the actual measured 
mis-alignments of the optics to produce a prediction of the 
final expected image quality. 
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